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ABSTRACT 

This  paper  investigates  the  feasibility  of  an  air-to- 
surface  missile  delivery  system  employing  impact  prediction 
to  assist  a  human  operator  in  controlling  the  system.   The 
impact  prediction  information  was  used  to  drive  a  display 
which  the  operator  used  in  controlling  the  flight  of  the 
weapon.   It  was  desired  to  find  whether  impact  orediction 
proved  useful  in  weapons  control. 
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I .   INTRODUCTION  AND  STATEMENT  OF  THE  PROBLEM 

A.   TYPICAL  ASM  WEAPONS  SYSTEMS 

Most  older  air-to-surface  weapons  share  the  characteris- 
tic that  once  the  weapon  is  launched  it  is  under  the  influence 
of  only  the  environment,  e.g.,  gravity,  wind  drift,  drag,  etc. 
Examples  in  this  category  of  weapons  are  bullets,  general- 
purpose  bombs,  and  rockets.   This  makes  it  the  responsibility 
of  the  weapons  officer  to  insure  that  the  fire-control  solu- 
tion is  accurate  enough  for  the  weapon  to  impact  at  the  tar- 
get within  the  kill  radius  of  the  warhead.   In  the  normal 
tactical  situation  this  is  a  difficult  task  to  perform  prop- 
erly.  As  a  result  some  of  the  newer  weapons  systems  have 
been  designed  with  a  capability  for  guidance  of  the  weapon 
to  the  target. 

These  newer  systems  require  that  to  provide  guidance  for 
the  weapon,  the  weapons  control  officer  must  maintain  visual 
contact  with  the  weapon  during  its  flight.   Again,  this  solu- 
tion to  the  problem  suffers  from  tactical  restrictions.   In 
a  combat  environment  it  is  very  difficult  to  maintain  a  flight 
path  which  allows  the  weapons  control  officer  to  maintain 
visual  contact  with  the  weapon  during  its  flight. 

In  order  to  overcome  the  shortcomings  of  these  systems, 
even  newer  systems  such  as  the  Phoenix  employ  very  sophisti- 
cated missiles  which  allow  the  launching  aircraft  to  stand-off 
from  the  target  and  still  deliver  a  weapon  so  as  to  destroy 
or  seriously  damage  the  target.   One  aspect  of  these  systems 


which  is  objectionable  is  the  high  cost  of  the  missiles  which 
are  employed.   As  a  result  this  study  was  made  of  a  system 
which  would  employ  a  stand-off  delivery  capability  while  em- 
ploying a  relatively  unsophisticated,  and  therefore  inexpen- 
sive missile  to  deliver  the  warhead. 

B.   ANALYSIS  OF  STAND-OFF  WEAPONS  DELIVERY 

Stand-off  weapons  delivery  is  best  defined  as  the  ability 
to  remain  outside  the  destructive  range  of  the  enemy's  weapons 
and  deliver  ordnance  on  the  enemy  installation.   The  advan- 
tages of  remaining  outside  the  destructive  envelope  of  the 
enemy's  weapons  are  obvious.   The  disadvantages  of  standing 
off  from  the  target  are  numerous.   The  greater  the  stand-off 
distance,  the  more  accurate  the  guidance  and  system  signals 
must  be  as  well  as  the  knowledge  of  the  target  location. 
This  is  necessary  to  provide  a  reasonable  probability  of  tar- 
get kill. 

In  its  most  basic  form  the  air-to-surface  weapons  delivery 
problem  consists  of  two  parts:   the  detection  and  location  of 
the  target  by  some  means,  and  delivery  of  a  weapon  to  the 
target. 

The  solution  which  is  proposed  in  this  paper  is  based  on 
three  assumptions  about  what  the  author  feels  an  air-to- 
surface  weapons  system  should  incorporate.   These  three  assump- 
tions are  1)  that  the  target  should  be  located  by  some  method 
which  does  not  endanger  the  launching  aircraft;   2)  that  the 
primary  purpose  of  air-to-surface  missiles  is  to  deliver  ord- 
nance on  the  target,  not  to  carry  a  large  amount  of  expensive 
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electronic  equipment  to  be  blown  to  pieces  along  with  the 
target;  and  3)  that  to  improve  the  overall  accuracy  of  the 
system,  the  missile  will  be  under  some  form  of  command  guid- 
ance during  its  flight  to  the  target. 

The  solution  which  is  proposed  in  this  paper  is  one  of 
the  "shoot  and  loiter"  school  of  weapons  control.   Those  of 
this  line  of  thought  contend  that  it  is  more  economical  to 
shoot  a  relatively  unintelligent  missile  at  the  target  and 
maintain  the  sophisticated  system  parts  in  the  launching  air- 
craft.  This  requires  that  the  launching  aircraft  remain  in 
the  target  area  until  the  weapon  has  hit  the  target  which 
increases  the  probability  of  the  launching  aircraft  being 
destroyed. 

The  opposite  point  of  view  is  that  the  aircraft  should 
launch  the  missile  and  get  out  of  the  target  area.   This 
approach  necessitates  the  use  of  a  very  sophisticated  missile 
but  reduces  the  probability  of  the  launching  vehicle  being 
destroyed. 

The  trade-off  to  be  considered  between  the  two  approaches 
is  to  obtain  the  proper  balance  of  economy  versus  probability 
of  survival  of  the  launching  aircraft. 


II .   THE  PROBLEM:  ANALYSIS  AND  SOLUTION 

In  this  chapter  the  two  divisions  of  the  problem  will  be 
discussed  in  detail.   The  first  is  the  problem  of  target  loca- 
tion.  The  second,  delivery  and  guidance  of  the  weapon,  will 
be  subdivided  into  two  sections.   The  first  section  will  be 
concerned  with  the  initial  fire-control  solution  for  the  mis- 
sile flight  path  and  the  second  with  the  guidance  of  the 
missile  to  the  target. 

A.   LOCATING  THE  TARGET 

The  task  of  locating  the  target  will  be  limited  to  the 
case  of  a  radiating  target  such  as  a  radar.   This  limitation 
was  made  in  order  to  propose  a  solution  in  which  the  launch- 
ing aircraft  both  locates  and  destroys  the  target.   The 
assumption  of  a  radiating  target  allows  the  use  of  a  direction- 
finding  method  to  locate  the  target. 

Several  methods  of  direction  finding  are  available.   These 
include  triangulation  of  radiation  given  off  by  the  target 
using  two  or  more  receivers  at  known  locations,  the  use  of 
filtered  sequential  measurements  of  bearing  and  depression 
angles  from  a  single  aircraft  at  known  positions,  or  even  the 
simple  method  of  having  two  ground  observers  locate  the  tar- 
get.  For  this  solution  the  second  method  was  used. 

1 .   Basic  Geometry  of  Target  Location 

In  Fig.  1  the  basic  principle  of  location  by  triangu- 
lation is  illustrated.   The  points  1  and  2  could  be  either 
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Figure  1.   Target  Location  by  Triangulation 
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Figure  2.   Definitions  of  Azimuth  and  Depression  Angles 
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two  distinct  receivers  listening  to  a  signal  at  the  same  time, 
or  they  could  be  the  two  locations  of  a  single  aircraft  at 
different  times. 

Figure  2  gives  a  definition  of  the  angles  involved  in 
all  the  target-location  schemes.   The  azimuth  angle  is  that 
angle  measured  from  a  true  north  reference  positive  clockwise 
to  tne  line-of-sight  (LOS)  between  the  aircraft  and  the  target, 
The  depression  angle  is  the  angle  measured  from  a  reference 
plane  normal  to  a  line  connecting  the  aircraft  and  the  center 
of  the  earth.   An  angle  which  rotates  below  the  normal  plane 
is  considered  negative  and  one  which  rotates  above  the  plane 
is  considered  positive.   Appendix  E  contains  a  discussion  of 
the  reference  systems  used  in  the  entire  problem. 

In  the  following  two  sections  two  different  methods 
of  location  of  the  target  by  use  of  angular  information  will 
be  presented. 

2 .   Location  Using  Positional  and  Angular  Information 

The  first  of  the  two  methods  for  target  location  will 
be  explored  here  for  the  two-dimensional  case  to  show  the 
principle  involved.   An  analysis  of  the  three-dimensional 
case  is  contained  in  Appendix  A. 

Figure  3  shows  the  geometry  of  the  two-position  method 
of  target  location.   The  aircraft,  at  two  different  times  and 
at  different  positions,  measures  the  azimuth  angle  from  the 
aircraft  to  the  target.   At  the  same  time  as  the  angle  is 
measured,  the  position  of  the  aircraft  is  recorded.   Assume 
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that  the  target  coordinates  are  X  and  Y , .   The  measurements 
are  listed  below. 

Time      Aircraft  Location      Azimuth  Angle 
1  X,  Yn  6, 


X. 


e. 


From  Figure  3  it  can  be  seen  that 

X  -X 
Tan  6   =  r^T^ 
1    *t  *1 


X  -x2 

Tan  9   =  — -T- 
xt   2 


(II-l) 


(H-2) 


Writing  the  two  equations  as  simultaneous  equations 
and  solving  them: 


X   -  X1   =    (Yt-Y,)  Tan  91 


Xt  -  X2  =  (Yt-Y2)  Tan  ^ 


Xt  -  Y   Tan  9,  =  X1  -  Y   Tan  61 


X   -  Y   Tan  e2  =  X2  ~  Y2  Tan  92 


(II-3) 


(II-4) 


Letting  A  =  Tan  9,  ,  B  =  Tan  9,,  C  =  X,  -  Y,  Tan  6,, 
and  D  =  X~  -  Y~  Tan  9~  and  writing  the  equations  in  matrix 
form 


1   -A 


1   -B 


X, 


(H-5) 


Using  normal  matrix  algebra  techniques,  the  solution 
can  be  shown  to  be 

(H-6) 


DA-BC 
Xt     A=B~ 


y  =  5z£ 

t    A-B 


(H-7) 
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Expanding  the  results  into  their  full  form  by  replacing 

the  A,  B,  C,  and  D  with  their  respective  values  gives 

TanB-j^  (X2-Y2Tan62)  -Tan62  (X.-Y.Tane, ) 
Xt  =  Ta"n6n-Taneo  (H-8) 


and 


(X  -Y9Tan69) - (X, -Y  Tan9, ) 
*t  "  TanB.-Tane^        *         K±x    *> 


This  technique  of  locating  the  target  can  be  quite 
effective  and  accurate  if  the  two  positions  used  in  deriving 
the  solution  are  far  enough  apart  and  if  the  bearing  angle 
measurements  are  respectably  accurate.   Ideally  the  two  bear- 
ing angles  should  be  separated  by  ninety  degrees.   This  sep- 
aration will  reduce  to  a  minimum  the  effects  of  any  noise  in 
the  measurement  of  the  angles.   If  the  target  radiates  for 
only  short  periods  of  time  and  the  ambient  conditions  are 
such  that  the  measurements  are  very  noisy,  the  derived  target 
location  could  be  greatly  in  error.   For  this  reason  an  al- 
ternate approach  to  the  problem  of  locating  the  target  was 
taken.   The  alternate  method  uses  only  angular  information 
and  a  hybrid  filtering  technique. 

3 .   Location  Using  Angular  Information  and  Hybrid  Filtering 

The  target  location  method  discussed  in  this  section 
is  the  result  of  research  done  by  Demetry  [1] .   The  method  em- 
ploys an  unusual  hybrid  filtering  technique  in  which  a  Kalman 
filter  is  considered  linear  for  calculation  of  a  gain  sched- 
ule but  in  which  prediction  is  done  nonlinearly.   The  hybrid 
filter  will  be  discussed  here  briefly. 
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It  is  desired  to  locate  a  radiating  target  by  a  method 
in  which  only  a  bearing  angle  to  the  target  is  measurable. 
The  equations  below  give  the  relationships  used  in  the  formu- 
lation of  the  filter. 

2klk  "  *(§k|k'§k|k'Ik>  (II-10) 


where 


a  = 


XT 
YT 


k|k 
klk 


(11-11) 


is  the  estimated  target  location  at  time  k  and 


k    - 


rx(t)i 

x(t) 
y(t) 
y(t) 


(n-12) 


is  the  aircraft  vector. 


The  aircraft  position  and  velocity  components  are 
considered  to  be  measurable  without  noise  for  this  evalua- 
tion as  they  are  throughout  the  remainder  of  the  problem. 

The  bearing  angle  dynamics  are  described  by  the 
differential  equations 


^(t)  = 


e(t) 

S(t) 

e(t) 

=  f (a,C(t) ,t) 


(11-13) 


where  f  is  a  nonlinear  vector  function.   The  angle  measurement 
at  any  discrete  time  k  is 


Zk  "  H^k  +  Vk'  k  -  T 


(11-14) 


in  which  H  =  [10  0]  and  v,  is  the  error  in  the  bearing  angle. 
The  v,  is  random  in  nature  with  an  assumed  mean  value  of  zero. 

The  gain  schedule  was  calculated  by  assuming  linear 
dynamics  and  constant  acceleration  for  the  process  and  using 
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the  recursive  Kalman  gain  equations  [2].   However,  the  predic- 
tion of  angular  information  was  done  by  using  the  actual  angu- 
lar dynamics  in  a  discretized  form  rather  than  employing  the 
linear  <t>  matrix  which  was  used  in  computing  the  filter  gains. 
These  equations,  by  using  true  system  dynamics  to  make  approx- 
imate predictions,  reduce  the  location  error  since,  in  effect, 
they  consider  the  system  linear.   In  this  system  this  is  not 
a  bad  assumption  since  the  sample  period  is  relatively  short 
compared  to  the  rates  of  change  of  the  system  variables. 


r 


§k+nk =  tan_1 


XTk i k-*k 
%ik-yk-ykTi 


(H-15) 


;  ^k[XTklk-xk] 

e       |     =  — — ^ — £- di-16) 

1         fYTklk-^kT]2+[XTk|k-xk]2 


2^k2[^klk-xk][^klk-^kT] 


(H-17) 


k+1'k        t[%|k-yk-ykT]^^k|k-xk]2}2      • 

This  is  a  brief  discussion  of  this  particular  method 
of  location  of  a  target.   For  the  reader  interested  in  pursu- 
ing this  further,  Ref.  1  will  be  of  interest. 

B.   THE  FIRE-CONTROL  SOLUTION 

Before  it  is  possible  to  launch  a  missile  at  a  target  it 
is  necessary  to  develop  a  solution  for  a  flight  path  which 
will  take  the  missile  from  the  launch  point  to  the  target. 
The  accuracy  with  which  the  solution  need  be  computed  is  en- 
tirely a  function  of  the  amount  of  guidance  which  will  be 
available  to  influence  the  missile  during  its  flight.   Also, 
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before  an  initial  solution  can  be  computed,  it  is  necessary 
to  know  what  approximate  flight  path  the  missile  will  follow 
from  the  launch  point  to  the  target. 

In  the  solution  which  is  being  investigated  here ,  the 
proposed  flight  path  for  the  missile  is  a  nominal  parabolic 
path.   The  path  the  missile  actually  travels  will  not  be  para- 
bolic due  to  air  drag,  control  effects,  and  other  outside 
influences . 

There  are  two  main  reasons  for  selecting  a  parabola  as  a 
nominal  flight  path.   The  first  concerns  a  relationship  be- 
tween the  height  of  a  parabola  and  distance  to  the  target 
from  a  point  under  the  peak  of  the  parabola.   This  ratio  is 
referred  to  as  the  A/R  ratio  and  its  importance  is  discussed 
in  Appendix  B.   The  second  consideration  is  that  the  missile 
will  be  under  the  influence  of  gravity  in  a  very  helpful  man- 
ner.  The  gravitational  force  will  help  to  maintain  the  veloc- 
ity of  the  missile  and  hence,  the  controllability  of  the 
missile,  after  the  boosted  phase  of  the  missile's  flight. 

1 .   Geometry  of  the  Fire  Control  Solution 

The  general  geometry  of  the  fire  control  problem  is 
illustrated  in  Fig.  4.   In  the  solution  the  aircraft,  point 
A/C  on  the  figure,  is  considered  to  be  at  the  X  origin  of  the 
figure  and  the  target,  point  T  in  the  figure,  is  considered  to 
lie  on  the  X-axis.   The  parabola  will  be  defined  such  that  it 
contains  the  two  points  of  the  launching  aircraft  (A/C)  and 

the  target  (T) .   Since  a  parabola  is  defined  by  a  second-order 

2 

polynomial,  e.g.,  of  the  form  y  =  ax  +  bx  +  c ,  some  third 

point  on  the  parabola  must  be  defined. 
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Figure  4.   The  Nominal  Parabolic  Flight  Path 
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The  Nominal  Parabolic  Flight  Path  with 
Distances  Defined 
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In  Fig.  5,  the  parabola  of  Fig.  4  is  redrawn,  this 
time  witn  three  points  and  four  distances  defined.   For  the 
purpose  of  defining  the  nominal  parabola,  the  launching  air- 
craft is  placed  on  the  Z  axis  so  that  range  is  zero  at  the 
origin.   The  altitude  of  the  aircraft  is  labeled  zA.   A  sec- 
ond point  is  the  location  of  the  target.   The  target  is  at 
zero  elevation  and  a  distance  D  from  the  origin.   The  third 
point,  which  is  necessary  for  the  parabola  to  be  defined,  is 
the  peak  of  the  parabola.   The  altitude  of  this  point  is  A 
and  it  is  D  -  R  from  the  origin.   Using  these  three  points 
it  is  possible  to  solve  for  the  coefficients  of  the  second- 
order  polynomial  defining  a  parabola. 

Writing  the  equations  of  the  parabola  at  the  three 
points  gives 

a(0)2  +  b(0)  +  c  =  zA  (11-18) 

a(D-R)2  +  b(D-R)  +  c  =  A  (11-19) 

a(D)2  +  b(D)  +  c  =  0  (11-20) 

R  =  kD    .  (11-21) 

Equation  11-21  is  used  to  define  the  relationship 
between  the  range  from  the  origin  to  the  target  and  the 
range  from  under  the  peak  of  the  parabola  to  the  target. 
2 .   Method  of  Obtaining  the  Launch  Angle 

In  a  system  in  which  the  launch  point  and  target  point 
are  coplanar  in  the  earth  plane,  a  ballistic  mass  will  travel 
furthest  when  the  velocity  vector  is  initially  directed  45 
degrees  above  the  horizontal.   This. is  true  as  long  as  the 
mass  has  no  external  forces  acting  on  it  other  than  gravity. 
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Two  factors  involved  in  this  investigation  make  the 
choice  of  a  45-degree  launch  angle  unacceptable.   The  first 
is  tnat  the  launching  point  and  target  are  not  earth  coplanar 
and,  secondly,  that  the  mass  is  subject  to  thrust  from  a 
rocket  motor.   Thus,  the  problem  of  computing  the  launch  angle 
becomes  one  of  maximizing  the  range  as  a  function  of  the  dis- 
tance between  the  target  and  launching  point  in  the  Z  plane. 

The  results  of  this  particular  maximization  were 
interesting.   For  aircraft  locations  which  were  ud  to  5000 
meters  noncoplanar  with  the  target,  the  launch  angle  for  max- 
imum distance  was  found  to  vary  only  one  degree,  between  4  7 
and  48  degrees.   (The  launch  angle  as  used  here  is  defined  as 
the  angle  between  the  velocity  vector  and  the  horizontal  at 
the  launch  point.   This  is  not  the  angle  of  the  velocity  vec- 
tor at  the  time  of  burnout  which  is  the  angle  which  deter- 
mines the  maximum  distance  of  flight.)   The  launch  angle  was 
set  at  a  constant  47.5  degrees  as  a  result.   This  angle  should 
give  the  maximum  range  regardless  of  range  to  the  target. 
Starting  with  the  launch  angle  that  gives  maximum  range  should 
prove  entirely  acceptable  since  one  of  the  assumptions  about 
the  system  was  that  the  missile  would  be  under  some  type  of 
guidance  during  its  flight. 

C.   GUIDANCE  BY  THE  HUMAN  OPERATOR 

Tnroughout  this  paper  mention  has  been  made  of  the  fact 
that  the  missile  will  be  guided  by  the  human  operator  through- 
out its  flight.   In  actuality  the  missile  was  not  released  to 
the  operator's  control  until  after  it  had  finished  the  boosted 


21 


phase  of  its  flight.   During  this  time  the  system  turned  the 
missile  up  to  the  launch  angle  and  yawed  it  to  the  proper 
flight  heading.   Once  the  missile  started  its  free  flight, 
a  display  of  the  predicted  impact  point  of  the  missile  rela- 
tive to  the  target  was  generated  for  the  missile  control 
officer  (MCO) .   At  the  same  time,  the  control  surfaces  of  the 
missile  were  set  to  a  neutral  position  and  control  of  the 
missile  was  given  to  the  MCO  who  then  guided  the  missile  the 
remainder  of  the  distance  to  the  target. 

In  order  for  the  predicted  missile  impact  point  to  be 
displayed,  it  was  first  necessary  to  track  the  missile  and 
based  on  the  tracking  information,  predict  its  impact  point. 
The  prediction  information  was  then  presented  to  the  MCO  in 
such  a  manner  that  he  could  make  decisions  as  to  the  guidance 
control  to  be  given  to  the  missile. 

1 .   Tracking  the  Missile  in  Flight 

In  this  problem  the  missile  was  tracked  in  flight  by 
a  radar  which  was  capable  of  providing  four  items  of  informa- 
tion.  These  were  the  range  to  the  missile,  the  range  rate  of 
change  to  the  missile,  the  azimuth  angle  to  the  missile,  and 
tne  depression  angle  to  the  missile.   The  reference  system  by 
whicn  the  angles  were  measured  is  discussed  in  Appendix  E. 
The  range  was  measured  in  the  positive  sense  while  the  range 
rate  was  considered  positive  for  a  missile  whose  range  was 
increasing  and  negative  for  a  missile  which  was  closing  the 
launching  aircraft. 
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Since  the  radar  was  operating  in  a  natural  environ- 
ment, the  information  provided  by  the  radar  was  not  perfect. 
This  noise  would  introduce  errors  in  the  estimation  of  the 
location  of  the  missile  durinq  its  flight,  so  it  was  desirable 
to  use  some  technique  to  reduce  the  effect  of  the  noise  on 
the  measurements.   The  technique  used  was  Kalman  filtering. 
Before  deriving  the  necessary  Kalman  filter  for  this  problem, 
the  desired  missile  information  vector  should  be  investigated. 

Figure  6  presents  the  geometry  for  developing  the 
missile  information  vector  from  the  spherical  vector  of  the 
radar.   The  desired  missile  information  vector  contains  six 
components:   the  XYZ  cartesian  location  of  the  missile  rela- 
tive to  the  aircraft  and  the  XYZ  cartesian  velocities,  again 
relative  to  the  aircraft.   The  acceptability  of  providing 
only  a  six-element  state  vector  for  the  missile  is  dependent 
upon  the  sample  rate  at  which  the  state  vector  is  developed. 
Even  though  the  missile  is  being  acted  upon  by  external  forces, 
it  is  possible  to  consider  the  missile  to  be  a  piece-wise 
constant-velocity  missile  in  either  spherical  or  cartesian 
coordinates  if  the  interval  between  samples  is  sufficiently 
small . 

From  Fig.  6  the  equations  for  the  missile  location 
relative  to  the  aircraft  can  be  developed.   Then  by  the  use 
of  the  time  derivatives  of  the  equations  of  location,  the 
equations  of  the  cartesian  velocities  can  be  found.   First, 
the  equations  of  location. 
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Z   =  R  sin  <J>  (11-22) 

m 

X   =  R'  sin  9  =  R  cos  <$>    sin  e  (11-23) 

m 

Y   =  R1  cos  6  =  R  cos  <J)  cos  6  (11-24) 

m 

For  the  case  shown  in  Fig.  6,  the  angle  tf>  is  negative, 
giving  the  proper  negative  value  for  Z  .   Likewise,  the  signs 
of  the  remaining  two  values  will  be  correct  depending  on  the 
quadrant  of  the  angle  0. 

To  form  the  equations  of  the  velocities  in  the  carte- 
sian frame,  it  is  necessary  only  to  take  the  time  derivatives 

of  equations  11-22,  11-23,  and  11-24.   These  are  given  below. 

•  •    • 

Z   =  R  cos  6    d>  -  R  sin  $  (11-25) 

m  r  r 

.  .  •      • 

X      =    R   cos    cj>    cos    0    0    -    R   sin    <b    sin    0    6   +    R   cos    d>    sin    0 
m 

(11-26) 
•  •  •  . 

Y      =    -R   cos    4>    sin    0    0    -    R   sin    cj>    cos    0    cj>   +    R   cos    <J>    cos    0 

(11-27) 
From  equations  11-22  through  11-27  a  list  of  the 
necessary  spherical  values  of  missile  movement  to  develop 
the  cartesian  state  vector  can  be  collected.   These  equations 
show  that  it  is  necessary  to  have  a  vector  of  spherical  values 
of  range,  azimuth,  depression,  and  the  associated  rates  of 
change  of  the  three  quantities.   From  this  requirement  a  Kal- 
man  filter  can  be  developed  which  will  not  only  reduce  the 
effect  of  the  measurement  noise  on  the  quantities  which  are 
being  measured,  but  also  provide  estimates  of  the  angular 
rates  which  are  not  measured. 

The  Kalman  filter  for  filtering  and  estimating  the 
radar  vector  has  the  general  Kalman  predictor-corrector  form 
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of  equation  11-28.   The  G,  is  computed  from  the  recursive 
equations  of  11-29,  11-30,  and  11-31.   More  complete  details 
of  the  derivation  of  the  filter  and  the  values  used  in  the 
gain  schedule  computation  are  given  in  Appendix  C. 

iik  ■  xjcik-i +  Gk<zk-H4ik-i'  (ii-28) 


Gk  ■  pk|k-iHt(Hpk|k-iHt+Rrl  (ii-29» 

pk|k  -  (i-GkH)pk|k-i  (ii-30> 

Pk+l|k  ■  *Pk|k*fc  +  2  (II-31» 

2 .   Predicting  the  Missile  Impact  Point 

In  predicting  the  missile  impact  point  it  is  desirable 
to  have  an  algorithm  which  is  fairly  accurate.   This  is  the 
case  since  it  is  the  point  at  which  the  missile  is  predicted 
to  strike  that  will  be  used  by  the  MCO  to  develop  his  guidance 
commands.   A  requisite  increase  in  accuracy  over  a  ballistic 
model  can  be  obtained  by  introducing  the  effects  of  drag  into 
the  equations  of  motion  for  the  missile.   For  the  prediction 
algorithm  developed  here,  the  drag  was  considered  to  be  oro- 
portional  to  the  velocity  of  the  missile. 

In  this  solution  the  drag  coefficient  of  the  missile 
is  taken  to  be  constant  from  the  time  of  the  prediction  until 
the  missile  strikes  the  target.   This,  of  course,  is  an  invalid 
assumption  but  one  which  is  necessary  in  order  to  develop  a 
closed-form  solution  for  the  prediction  equations.   The  error 
which  is  introduced  by  this  assumption  will  have  the  effect 
of  causing  the  MCO  to  hold  the  missile  aloft  longer  than  would 
normally  be  necessary.   This  action  by  the  MCO  will  tend  to 
increase  the  A/R  ratio,  thus  making  the  missile  impact  on  the 
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target  more  vertically  and  hence  less  subject  to  error  from 
target  altitude  location  error.   (See  Appendix  B) . 

a.   Equations  of  Motion  for  Impact  Prediction 

The  equations  of  motion  for  impact  prediction  were 
written  for  each  of  the  three  component  directions  of  the 
cartesian  reference  system.   The  equations  were  based  on  an 
assumed  point-mass  model  of  the  missile  subjected  only  to 
drag  and  gravity  forces.   The  X  and  Y  equations  are  similar 
to  each  other,  varying  only  in  the  drag  coefficient.   The  Z 
equation  takes  into  account  the  effect  of  gravity.   The  equa- 
tions of  motion  are 

mX  =  -K  X  (11-32) 

my  =  -K2Y  (II-33) 

mZ  =  -K  Z  -  mG  (11-34) 


where 


Kx  =  p|x|Ar  (H-35) 

K0  *  p I Y I  A  (11-36) 

2    K  '  '  r 


m 


K3  *  p|Z|Ar  (II-37) 

and  m  =  missile  mass,  G  *  gravity,  and  A  =  reference  area. 

The  method  for  predicting  the  impact  point  of  the 
issile  is  to  first  compute  the  drag  coefficients  (K..  ,  K  , 
and  K-.)  which  are  then  considered  constant  until  impact.   The 
Z  motion  equation  is  then  solved  for  the  time  when  the  missile 
Z  location  is  equal  to  the  target  Z.   First,  it  is  necessary 
to  solve  equation  11-34  for  the  missile  Z  as  a  function  of 
time . 
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b.   Impact  Prediction  Equations 

The  result  of  solving  for  the  Z  motion  as  a  func- 
tion of  time  is 


Z(t)  =  (1-A)e"at  -  Bt  +  C  (11-38) 


where 


A=  4^-  +  ^-  (H-39) 

K3      K7 

B  =  p-  (11-40) 

K3 

C  =  Z  (0)  (11-41) 

m 

K3 
a   =  —   .  (11-42) 

m 

This  mixed  equation  is  not  directly  solvable  for 

that  value  of  time  such  that  Z   is  equal  to  Z,  .   The  Newton- 

m     M        t 

Raphson  iterative  method  is  employed  to  find  the  value  of 

time  which  will  satisfy  the  equation  for  a  specific  set  of 

conditions  (a  particular  Z  (0)  and  Z  (0) ) .   The  value  of  time 

mm 

thus  found  is  substituted  into  the  X  and  Y  solutions  to  find 
the  predicted  impact  point  of  the  missile.   The  X  and  Y 
predicted  positions  at  impact  are  given  by  equations  11-43 

and  11-44. 

mX  (0)      -K.t./ni 

X^  =  X  (0)  +  -JS (l-e  l    1   )  (11-43) 

p    m       K.. 

mY  (0)      -K„t./m 

Y^  =  Y  (0)  +  -t=2 (l-e  Z    X   )  (11-44) 

p    m       K» 

where  t.  is  the  value  of  impact  time  from  the  solution  of 
equation  11-38. 
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3.   Display  of  the  Situation 

In  order  for  the  MCO  to  make  intelligent  decisions 
regarding  the  application  of  control  to  the  missile  during 
flight,  there  had  to  be  some  method  for  displaying  the  pre- 
dicted results  of  the  current  flight  path.   What  was  displayed 
to  the  MCO  was  the  current  prediction  of  the  impact  point 
of  the  missile  assuming  that  there  would  be  no  further  control 
exercised  over  the  missile.   The  display  is  generated  such 
that  the  center  of  the  screen  corresponds  to  the  last  esti- 
mate of  the  target  position  and  the  predicted  impact  point  is 
displayed  relative  to  the  center  of  the  screen.   The  disolay 
is  also  rotated  so  that  motion  of  the  impact  ooint  will 
correspond  to  control  stick  motion.   Figure  7  illustrates 
the  display  screen  and  Appendix  H  further  discusses  the  man/ 
system  interface. 
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III.   THE  SIMULATION  OF  THE  PROBLEM 

A  proposed  solution  to  a  problem  cannot  be  considered  to 
be  valid  unless  it  is  tested  in  some  way  in  practice.   One 
method  of  evaluation  is  to  build  an  actual  system  and  test 
it  in  a  practical  operating  environment.   This  is  generally 
too  expensive  to  be  feasible.   Another  method  of  testing  a 
proposal  is  to  simulate  the  system  on  an  appropriate  computer, 
The  second  method  was  chosen  to  evaluate  the  problem  solution 
proposed  in  this  paper. 

This  simulation  was  done  on  two  different  computers.   The 
first  was  an  IBM  System  360  Model  67,  operating  under  the  MVT 
operating  system.   This  system  was  used  for  initial  program 
debugging  since  the  Operating  System  has  an  excellent  error- 
monitoring  system.   The  second  system  was  a  hybrid  computer 
consisting  of  an  SDS  9300  digital  computer  and  a  COMCOR  CI- 
500  0  analog  computer. 

A.   SIMULATION  OF  TARGET  LOCATION 

The  simulation  of  the  target  location  was  accomplished 
by  using  the  results  of  work  in  passive  target  location  done 
by  Demetry  [1] .   The  results  of  the  hybrid  location  method 
discussed  in  Chapter  II  were  used  by  making  use  of  the  calcu- 
lated mean-square  error  values  in  target  X  and  Y  location 
as  a  function  of  time,  from  the  hybrid  location  algorithm. 

Using  these  values  involved  a  process  of  both  normaliza- 
tion and  interpolation.   Since  the  simulation  is  designed  for 
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a  time  increment  of  one  second,  the  target  location  mean- 
square  errors  had  to  be  interpolated  since  they  were  computed 
on  a  time  increment  of  two  seconds  between  observations. 

Tue  interpolation  was  made  on  the  assumption  that  the 
mean-square  errors  of  location  were  piece-wise  linear  over 
tne  measurement  interval  involved.   The  normalization  process 
was  necessary  since  when  the  values  were  generated,  the  air- 
craft was  flying  a  set  course  which  had  a  direct  effect  on 
tne  value  of  the  mean-square  error.   Thus,  the  values  were 
normalized  by  multiplying  them  by  a  factor  which  was  the 
closest  point  of  approach  (CPA)  to  the  target  for  the  prob- 
lem divided  by  the  CPA  to  the  target  in  the  mean-square-error 
computation  program. 

B.       SIMULATION    OF   THE    FIRE-CONTROL    SOLUTION 

A  real-time  simulation  of  the  fire-control  solution  was 
eliminated  after  some  research  into  the  problem.   A  constant 
launch  angle  was  used  in  its  place.   The  solution  was  arrived 
at  oy    assuming  the  missile  to  be  a  drag-free  ballistic  mass 
with  the  propulsion  system  of  the  true  missile  attached. 
Tnis  approximation  was  used  in  an  iterative  method  of  comput- 
ing the  angle  which  would  give  the  maximum  flight  down  range 
for  a  particular  launch  altitude.   The  iteration  was  performed 
on  both  the  launch  altitude  and  the  launch  angle.   For  this 
investigation  the  angular  increment  was  taken  as  one  degree 
and  the  altitude  increment  was  fifty  meters. 

For  each  launch  altitude,  the  angle  which  gave  the  maximum 
down-range  flight  of  the  assumed  missile  was  plotted  against 
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the  launch  altitude.   The  results  were  surprising  in  that  it 
was  found  that  for  altitudes  between  fifty  and  five  thousand 
meters,  the  angle  varied  only  between  4  7  and  4  8  degrees.   The 
launcn  angle  is  defined  as  the  initial  angle  between  the  hor- 
izontal and  the  velocity  vector.   This  initial  angle  is  re- 
duced by  the  effect  of  gravity  to  an  angle  at  burnout  which 
will  give  the  maximum  down-range  flight  of  the  mass.   The 
force  of  gravity  is  thus  responsible  for  the  initial  angle 
being  greater  than  45  degrees. 

It  must  be  noted  that  this  angle  was  chosen  to  give  the 
maximum  down-range  flight  of  the  assumed  missile.   Thus,  by 
shooting  the  true  missile  at  this  angle  the  maximum  possible 
down-range  flight  would  be  possible  for  zero  control.   By 
the  application  of  control  to  the  missile  in  the  pitch  mode, 
this  range  can  be  either  increased  or  decreased.   In  the 
desired  situation  this  range  should  have  to  be  decreased  by 
using  down  pitch.   This  action  would  cause  an  increase  in 
tne  A/R  ratio  which  has  already  been  demonstrated  to  be 
desirable . 

C.   SIMULATION  OF  THE  MISSILE 

The  missile  model  was  the  most  important  part  of  the  simu- 
lation and  thus  was  treated  with  the  most  rigor.   The  attempt 
was  made  to  use  a  fairly  accurate  model  in  that  dynamics 
approximating  an  actual  missile  were  assumed.   Appendix  F 
contains  a  description  of  the  missile  model. 
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D.   GUIDANCE  SYNTHESIS 

The  guidance  synthesis  for  the  simulation  was  done  in  two 
parts:   one  part  for  the  simulation  while  it  was  executed  on 
tne  System  360  and  the  second  for  the  actual  data-gathering 
simulation  on  the  SDS/COMCOR  hybrid. 

1 .  System  360  Guidance  Synthesis 

The  guidance  synthesis  used  during  the  preliminary 
runs  on  tne  System  360  was  a  primitive,  bang-bang  control. 
Yaw  was  either  plus  or  minus  two  degrees  per  second  and  pitch 
was  a  plus  or  minus  0.5  radian  deflection  of  the  control  sur- 
face.  The  assumptions  on  which  this  control  was  synthesized 
were  tne  same  as  those  discussed  in  the  following  sections. 

2 .  SDS/COMCOR  Guidance  Synthesis 

Before  going  into  the  discussion  of  the  methods  used 
to  develop  the  guidance  commands  in  the  hybrid  simulation,  the 
assumptions  on  which  the  command  synthesis  was  based  will  be 
discussed. 

a.   Roll-Independent  Control 

The  assumption  that  control  is  roll  independent 
is  valid  since  the  majority  of  missile  control  systems  operate 
in  this  manner.   Also  the  assumption  allows  a  simplification 
of  the  command  vector  synthesis.   The  missile  is  thus  consid- 
ered to  either  be  roll  stabilized  or  the  control  system  in 
the  missile  is  so  designed  to  compensate  for  any  roll  of  the 
missile. 
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b.  Two-Component  Control  Vector 

With  roll  independence  assumed,  the  control  vector 
consists  of  two  components,  a  pitch  and  a  yaw  component. 
These  components  were  limited  in  their  magnitude  by  the  syn- 
thesis program. 

The  yaw  was  limited  to  a  maximum  rate  of  twenty 
degrees  per  second  by  the  control  mechanism  itself.   The 
pitch  was  limited  to  plus  or  minus  forty-five  degrees  of 
control  deflection.   Additional  limits  on  the  control  deflec- 
tions were  built  into  the  missile  model.   Yaw  rate  was  lim- 
ited so  that  radial  acceleration  on  the  missile  would  not 
exceed  ten  G's. 

c.  Control  Vector  Formation 

The  control  vector  was  formed  by  the  human  opera- 
tor directing  a  control  stick.   The  maximum  deflection  of  the 
stick  corresponds  to  the  maximum  amount  of  control  which 
could  be  generated.   This  control  was  then  applied  to  the 
missile  each  second  until  the  control  stick  was  neutralized. 
The  control  affected  the  missile  in  that  the  control  surfaces 
of  the  missile  remained  deflected  only  as  long  as  the  stick 
was  held  in  a  non-neutral  position. 
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IV.   DISCUSSION 

The  simulation  of  this  problem  was  performed  only  a  few 
times  with  the  completed  program.   As  a  result  it  was  not 
possible  to  draw  concrete  conclusions  about  the  validity  of 
tne  total  concept  tested.   However,  several  observations  were 
made  about  the  program  and  some  areas  where  further  investi- 
gation could  be  made  were  discovered. 

A.   OBSERVATIONS 

One  observation  from  the  running  of  the  program  supported 
tne  supposition  that  the  use  of  impact  prediction  in  an  air- 
to-surface  environment  was  effective.   The  errors  in  impact 
prediction  seemed  to  be  related  closely  to  errors  in  producing 
the  measured  missile  state  vector.   Overall  the  impact  pre- 
diction approach  to  missile  control  appeared  very  feasible. 

Another  observation  which  was  made  concerned  the  initial 
launch  angle  assumed  for  the  simulation.   Apparently  the  use 
of  a  fixed  angle  of  launch  intended  to  give  maximum  down- 
range  flight  was  not  as  good  an  assumption  as  first  thought. 
For  targets  which  were  close  to  the  launching  aircraft,  it 
was  possible  to  overfly  the  target  by  shooting  for  maximum 
range.   It  would  be  much  more  desirable  to  have  the  launch 
angle  be  a  function  of  the  initial  range  to  the  target. 

The  overfly  or  fly-by  effect  can  be  a  serious  problem. 
If  the  missile  bypasses  the  target  with  a  small  angular  offset 
to  either  side,  it  is  possible  that  immediate  application  of 
yaw  control  in  such  circumstances  will  cause  the  missile  to 
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circle  the  target  at  a  constant  range.   This  can  occur  since 
the  missile  is  limited  in  the  amount  of  radial  acceleration 
it  can  withstand.   It  should  be  possible  to  overcome  the  prob- 
lem by  the  use  of  tactical  doctrine  if  the  missile  does  over- 
fly the  target.   One  suggested  doctrine  would  be  to  continue 
the  fly-by  course  until  the  missile  is  far  enough  out  to  make 
the  turn  and  hit  the  target  area.   This  maneuver  is  illus- 
trated in  Fig.  8. 

B.   FURTHER  INVESTIGATION 

Three  areas  exist  which  the  author  feels  merit  further 
investigation.   The  areas  are  1)  prediction  accuracy  and 
methods;  2)  radar  filter  accuracy;  and  3)  display  methods. 

The  accuracy  of  impact  prediction  should  be  investigated 
further.   Although  it  was  shown  that  the  impact  prediction 
method  worked,  no  study  of  the  accuracy  of  the  method  em- 
ployed was  done.   Also,  it  is  possible  that  different  meth- 
ods of  impact  prediction  would  prove  to  be  more  accurate 
than  the  velocity-proportional  model  used  in  the  simulation. 

Prediction  accuracy,  from  all  indications,  was  very 
dependent  on  the  accuracy  with  which  the  missile  state  was 
estimated.   The  key  portion  of  the  program  in  the  missile- 
state  estimation  was  the  radar  vector  filter.   The  filter 
used  in  the  problem  was  based  on  the  assumption  of  a  constant- 
velocity  model  of  the  missile.   The  use  of  a  filter  based  on 
constant  acceleration  should  be  investigated  as  lags  were 
noted  in  the  measured  missile  state  compared  to  the  true  mis- 
sile vector  when  the  missile  was  making  hard  or  continued 

maneuvers . 
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The  method  of  displaying  the  prediction  information  to 
the  MCO  also  deserves  further  investigation.   The  display  as 
it  was  implemented  did  not  provide  any  indication  of  the  veloc- 
ity of  the  missile.   It  is  felt  that  if  the  MCO  were  given 
missile  velocity  information,  he  would  be  able  to  better  con- 
trol the  flight  of  the  missile.   This  is  especially  true  in 
the  fly-by  case. 

As  the  missile  bypasses  the  target,  the  display  gives  a 
false  indication  as  to  the  control  necessary  to  hit  the  tar- 
get.  This  problem  is  generated  as  a  result  of  the  almost 
180-degree  change  in  angle  by  which  the  display  is  rotated. 

Figure  9  illustrates  an  improved  display  which  should 
solve  the  rotation  problems  caused  by  fly-by  and  also  give 
velocity  vector  information  to  the  MCO.   The  display  is  of 
the  entire  area  of  the  problem.   The  axes  of  the  new  display 
are  neither  translated  nor  rotated  as  they  were  in  the 
simulation. 

Three  points  are  included  on  the  display.   One  is  the 
current  XY  measured  position  of  the  missile.   This  is  the 
left-most  point  in  Fig.  9.   The  other  two  points  are  the 
current  predicted  impact  point  and  the  last  estimate  of  the 
target  location.   A  line  from  the  missile  position  to  the  im- 
pact point  will  indicate  the  velocity  information  desired. 
Also  the  non-rotation  of  the  display  axes  should  alleviate 
the  false  indication  of  necessary  command. 
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Missile  Path 


Figure  8.   Course  for  Fly-by  Error  Correction 
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Figure  9.   Suggested  Improved  Display 
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C.   RESULTS  OF  A  RUN 

The  figures  in  this  section  show  the  results  of  only  one 
run  of  the  simulation  program.   The  figures  give  two  planar 
plots  of  the  flight  path.   One  is  a  plot  of  the  XY  plane 
flight  of  the  missile  (Fig.  10)  and  the  second  is  a  vertical 
plane  view  of  the  flight  path  (Fig.  11) . 

Other  plots  are  included  which  show  the  control  vector 
as  a  function  of  time  (Fig.  12)  and  the  missile  location 
errors  as  a  function  of  time  (Fig.  13-14) . 

It  should  be  emphasized  again  that  this  is  only  one  run 
and  should  not  be  considered  to  be  in  any  way  conclusive.   To 
reach  any  definite  conclusions  a  much  more  exhaustive  investi- 
gation would  have  to  be  made. 
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APPENDIX  A 

TARGET  LOCATION  BY  SUCCESSIVE  ANGULAR  MEASUREMENTS 
FOR  THE  THREE-DIMENSIONAL  CASE 

In  chapter  II  there  was  given  a  brief  description  of 
locating  a  target  in  the  two-dimensional  plane  by  processing 
two  sets  of  angular  measurements  and  two  associated  aircraft 
positions.   In  this  appendix  the  two-dimensional  case  will 
be  expanded  to  the  third  dimension. 

Figure  15  shows  the  three-dimensional  geometry  involved 
in  this  problem.   This  figure  is  then  shown  in  two  plane  views 
the  XY  plane  projection  in  Fig.  16  and  the  RZ  plane  projec- 
tion of  Fig.  17. 

From  Fig.  16  the  two  locations  of  the  launching  aircraft 
are  marked  as  X, ,  Y,  and  X~ ,  Y„  with  the  associated  azimuth 
angles  6,  and  9~.   Using  the  geometric  relationship  for  tan- 
gent, the  tangents  of  the  angles  shown  are 

X  -X 

v   v   =  tan  6  (A-l) 

t   1 

Xt~X2 

v   v   =  tan  9    .  (A-2) 

*t~x2         z 
Treating  the  equations  as  simultaneous  equations  and  setting 
them  up  for  solution, 

Xt  -  Y   tan  91  =  X1    -  Y^^  tan  61  (A-3) 

Xt  -  Y   tan  92  =  X2  -  Y2  tan  92  .        (A-4) 
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Letting  tan  9,  =  A,  tan  9   =  B ,  X   -  Y..  tan  9..  =  C,  and 

X2  -  Y 

tions  A-3  and  A-4  results  in 


^  „  tan  9?  =  D  and  substituting  these  values  into  equ 


a- 


Xt  -  AYt  =  C 


Xt  -  BYt  =  D  . 


(A-5) 
(A-6) 


Rewriting  these  equations  in  matrix  form  gives 

I   -A]    jit  1     -C- 


Ffc 


1   -B 

which  has  a  solution  of  the  form 

X,  " 


D 


(A-7) 


1   -A'  1  "C 


1   -B 


D 


(A-8) 


The  solution  of  this  equation,  arrived  at  by  normal  matrix 
algebra,  is 


(A-9) 


Making  substitutions  back  into  the  solution,  the  X,  and  Y. 
of  the  target  are 

tan  9   (X9-Y9tan99)  -  tan  99  (X-.-Y.tan6,) 

X^  =  ± ~ = r tt =-   (A-10) 

t  tan  9,  -  tan  9^ 


X^ 

-B 

A 

rC] 

t 

i 

A-B 

k 

-1 

1 

D 

_ 

_ 

. 

u          J 

Yt  - 


(X2-Y2tan92)  -  (Xj-Y^ane^ 

tan  9,  -  tan  9„ 


(A-ll) 


The  solution  for  locating  the  Z  coordinate  of  the  target 
will  now  be  derived. 

Figure  17  will  be  used  to  illustrate  the  derivation  of 
the  target  Z  location  problem.   Again  using  the  tangent 
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+Y 


X2(X2,Y2) 


X1(X1/Y1) 


+X 


Figure  16.   XY  Plane  View  Location  Geometry 


+Z 


+R 


Figure  17.   RZ  Plane  View  Location  Geometry 
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relationships  for  the  angles  involved,  it  can  be  seen  that 

Vzt 

-^ —  =  tan  (J)1  (A-12) 

Z2"Zt 

n    =  tan  <K   .  (A-13) 

U2 
From  Fig.  16  the  values  for  D,  and  D„  can  be  seen  to  be 

D1  =  ((Xt-Xx)2  -  (Yt-Y1)2)%  (A-14) 


2  2  ^ 

.        -  (Y  -Y 

2    v  v  t  2'  v  t   2 


D0  =  ((X  -Xor  -  (Y^-Yn)")2  .  (A-15) 


The  fact  is  very  evident  from  equations  A-14  and  A-15 
that  once  the  X  and  Y  coordinates  of  the  target  have  been 
found,  any  single  measurement  of  the  depression  angle  may 
be  used  along  with  the  Z  location  of  the  aircraft  at  the 
time  of  the  depression  angle  measurement  to  determine  the  Z 
coordinate  of  the  target. 
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APPENDIX  B 

EFFECT  OF  A/R  RATIO  AND  TARGET  Z  LOCATION  ERROR 
ON  RANGE  MISS  DISTANCE 

Several  references  have  been  made  in  this  paper  to  the 
A/R  ratio  and  its  importance  in  determining  missile  miss  dis- 
tances at  the  target.   This  appendix  is  devoted  to  a  discus- 
sion of  the  importance  of  this  effect. 

A.   DEFINITION  OF  A/R  RATIO  AND  ITS  EFFECT  ON  MISS  DISTANCE 

The  A/R  ratio,  as  defined  for  a  parabola,  is  the  ratio  of 
the  height  of  the  parabola  at  its  apex  to  the  distance  from 
a  point  under  the  apex  to  one  of  the  terminal  points  of  the 
parabola.   This  ratio  is  illustrated  in  Fig.  18.   It  can  be 
seen  from  the  figure  that  the  larger  the  A/R  ratio,  the  more 
nearly  perpendicular  the  terminus  of  the  parabola  at  its 
intersection  with  the  axis.   The  more  perpendicular  the  ter- 
minus of  the  parabola,  the  less  the  point  of  intersection 
with  the  target  plane  will  be  affected  by  an  error  in  locating 
the  target  plane.   Figures  19  and  20  illustrate  the  relation- 
ship just  described. 

In  Fig.  19,  which  illustrates  a  small  A/R  ratio,  it  can 
be  seen  that  for  an  error  in  target  Z  location,  AZ ,  a  rela- 
tively large  error  in  range,  AR,  is  generated.   The  alternate 
and  desirable  case  of  the  large  A/R  ratio  is  shown  in  Fig.  20. 
For  the  same  error  in  target  Z  location  as  in  Fig.  19,  the 
range  error,  AR,  is  much  smaller  than  in  the  small  A/R  ratio 
example . 
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Figure  19 .   Small  A/R  Ratio 


Figure  20.   Large  A/R  Ratio 
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B.   CURVES  RELATING  MISS  DISTANCE  AND  A/R  RATIO  AS  A  FUNCTION 
OF  TARGET  AZ. 

The  computation  of  a  family  of  curves  relating  miss  dis- 
tances and  A/R  ratio  as  a  function  of  target  AZ  is  straight- 
forward.  Assuming  a  ballistic,  drag-free  mass,  the  time  for 
that  mass  to  fall  the  distance  A  under  the  influence  of 
gravity  was  determined.   This  was  done  for  the  initial  condi- 
tion of  the  mass  having  a  zero  Z  velocity.   In  equation  form 

tk  =  ~-^  (B-l) 

K  Z 

2 

where  Z  =  -9.8  m/sec   and  A  is  the  initial  Z  of  the  mass  with 

the  target  Z  assumed  to  be  zero.  By  using  this  time  and  the 
knowledge  of  the  distance  R  from  the  apex  of  the  parabola  to 
the  target,  the  initial  R  velocity  was  computed  as 

I  Vr  =  b      ■  (B-2) 

1 

This  velocity  was  then  set  as  constant  throughout  the 

current  run  in  which  the  Z  location  of  the  target  was  varied 

away  from  the  perfect  zero  location  for  which  the  R  velocity 

was  computed.   For  each  value  of  AZ ,  a  time  to  fall,  tf,  was 

computed  and  then  substituted  into  the  equation 

R  =  V  t,   '.  (B-3) 

r  f 

The  difference  between  this  distance  thus  computed  and  the  R 
to  the  target  was  the  miss  distance.   The  series  of  computa- 
tions was  repeated  for  a  number  of  different  A/R  ratios  and 
AZ's. 

The  results  were  then  plotted  with  A/R  ratio  on  the  X 
scale  and  miss  distance  on  the  Y  scale.   The  family  of  curves 
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of  Fig.  21  was  generated  by  connecting  the  miss  distances  at 
the  different  A/R  ratios  which  were  the  result  of  the  same 
AZ. 

The  values  used  in  this  problem  were 

R  =  5000  meters   (constant) 

A's  such  that  A/R  ratios  of  .1,  .2,  .3,  .4, 

.5,  .6,  .7,  .8,  .9,  1.0,  1.5,  2.0,  2.5,  3.0, 

3.5,  and  4.0  were  generated. 

AZ's  used  were  1,  2,  5,  10,  20,  50,  100, 

200,  and  500  meters. 
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APPENDIX  C 

DERIVATION  OF  THE  KALMAN  FILTER 
FOR  THE  MISSILE-TRACKING  RADAR [2] 

In  chapter  II  the  Kalman  filter  for  the  system  radar  was 
discussed.   This  appendix  gives  the  assumptions  which  were 
used  in  the  derivation  of  the  filter. 

The  filter  was  designed  with  the  assumption  that  the 
missile-tracking  radar  provides  measurements  of  range,  range 
rate,  azimuth  angle,  and  depression  angle  to  the  missile 
relative  to  the  aircraft.   Another  assumption  was  that  the 
measurements  were  contaminated  by  additive  noise  whose  sta- 
tistics were  known.   It  was  also  assumed  that  the  missile  was 
not  subjected  to  any  random  perturbations.   Based  on  these 
assumptions,  some  of  the  matrices  of  equations  C-l,  C-3,  and 
C-5  may  be  defined. 

S       -  pkik-iHt(HPkik-iHt+Rrl  (c-1) 


pk|k    =  (I-GkH)pk|k-i  (c-2) 

Pk+l|k   =    ^klk*'  +  2  ((>3> 

*kik    =  4ik-i +  Gk(Zk-H*kik-i>  (c-4) 


^kik-i  =  *ik-iik-i  (c-5)1 


The  Q  matrices  are  identically  zero  since  it  was  assumed 
that  the  missile  was  subject  to  no  random  perturbations.   The 


The  symbols  used  in  the  equations  C-l  through  C-5  and 
elsewhere  in  this  appendix  differ  somewhat  from  the  svmbol- 
ogy  used  in  R.  C.  K.  Lee.  [Ref.  2]. 
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$    =0    =  $ 
rg     az     dp 


(C-6) 


R  matrices  are  not  zero  since  they  convey  to  the  filter  infor- 
mation regarding  the  variance  of  the  additive  noise  in  the 
measured  values. 

The  $  matrices  are  the  same  for  all  the  filters  since  the 
same  dynamics  (constant  velocity)  were  assumed  for  all  the 
filters.   Equation  C-6  gives  these  $  matrices  as 

'1   T" 

The  H  matrices  in  equations  C-l ,  C-2 ,  and  C-4  are  measure- 
ment matrices  which  determine  which  variables  are  measurable 
by  the  filters.   In  the  range  filter  H    is  a  2  X  2  identity 

matrix.   The  H    and  H,   are  1X2  matrices  with  the  (1,1) 

az       dp 

elements  being  unity. 

It  should  be  noted  that  equations  C-l,  C-2,  and  C-3  are 
data-independent  as  long  as  the  matrices  involved  in  computing 
them  are  data-independent  and  if  the  H,  R,  $ ,  and  Q  matrices 
are  constants.   If  this  is  true,  then  the  gain  schedule  can 
be  computed  offline  from  the  actual  simulation  of  the  problem 
or  operational  environment  and  simply  inserted  into  the  com- 
puter memory  for  use  at  the  proper  time.   With  the  precomputed 
gain  schedule  in  use,  the  implementation  of  the  filter  to 
process  data  in  real-time  is  a  relatively  easy  matter  and 
does  not  contribute  greatly  to  the  computational  load  on  the 
computer. 

An  important  consideration  in  designing  any  filter  is  the 
expected  value  of  the  measured  parameters  which  will  be  used 
to  initialize  the  filter.   This  is  especially  important  when 
choosing  values  for  the  initial  covariance  matrices,  P, iQ. 
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As  these  filters  were  to  be  employed  in  conjunction  with  a 
missile-tracking  radar,  the  assumptions  were  made  that  the 
initial  range  to  the  missile  and  range  rate  of  the  missile 
would  be  unknown  and  that  the  initial  angles  to  the  missile 
would  also  be  unknown.   These  assumptions  were  reflected  in 
the  covariance  matrix  initializations  when  the  gain  schedule 
was  computed.   The  values  used  are  given  for  the  range  and 
range  rate  filter  in  equation  C-7  and  the  angular  filters  in 
equation  C-8. 


10 


10 


10000.  m' 
0 

,2 


0 
1000.  m2/sec2 


0 


(C-7) 


2     2 

10  rad  /sec 


(C-8) 


10  rad- 
0 

The  values  used  for  the  R  matrices  in  the  gain  schedule 
computation  were 

n100  m2       0 
0 


rg 


10  m  /sec 


(C-9) 


R   =  R ,   =  .0004  rad' 
az     dp 


(C-10) 


By  using  very  large  initial  covariance  matrices  and  set- 
ting the  initial  estimates  of  the  state  vectors  to  zero,  the 
gain  element  (1,1)  for  the  first  time  increment  will  be  very 
close  to  unity,  thus  causing  the  first  measurement  of  the 
particular  state  variable  to  become  the  first  estimate  of 
that  variable.   As  time  passes,  however,  the  incoming  meas- 
urement is  weighted  less  and  less  as  the  gain  goes  down,  i.e., 
the  filter  becomes  more  confident  in  the  accuracy  of  its 
estimates . 


59 


However,  this  reduction  of  gain  as  a  function  of  time  can 
sometimes  be  self-defeating.   In  the  case  in  which  this  filter 
was  employed,  the  reduction  of  the  gain  with  time  could  be 
disastrous.   Some  corrective  action  had  to  be  taken  to  prevent 
the  gains  from  approaching  zero  and,  as  a  result,  having  the 
filter  ignore  new  data.   Otherwise,  as  the  missile  maneuvered, 
the  filters  would  have  very  quickly  begun  to  produce  very  bad 
estimates  of  the  radar  state  vector. 

Thus,  the  effect  which  is  desired  is  one  which  will  not 
allow  the  gains  to  approach  zero.   There  are  two  actions  which 
will  prevent  this  from  occurring:   letting  the  0  matrices  be 
non-zero,  or  truncating  the  gain  schedule  after  some  length 
of  time  kT .   Since  the  Q  matrices  have  already  been  considered 
to  be  zero,  there  remains  gain-schedule  truncation.   By  trun- 
cating the  gain  schedule  the  filters  continue  to  accept  in- 
coming measurements  as  containing  significant  information, 
and  the  amount  of  core  storage  in  the  computer  is  reduced 
greatly . 

The  gain  schedules  in  this  particular  simulation  were 
truncated  in  such  a  manner  that  the  previous  estimate  produced 
by  the  filters  were  weighted  by  a  factor  of  approximately  0.8 
and  the  incoming  measurements  were  weighted  by  a  factor  of 
approximately  0.2. 
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APPENDIX  D 

DERIVATION  OF  THE  SOLUTION  TO  THE 
IMPACT  PREDICTION  EQUATIONS 

This  appendix  discusses  the  derivation  of  the  impact 
prediction  equations  used  for  the  missile  impact  point  pre- 
diction.  Two  methods  are  discussed:   one  in  which  the  drag 
force  is  proportional  to  the  velocity  and  a  second  in  which 
the  drag  force  is  proportional  to  the  velocity  squared.   The 
first  was  used  in  the  simulation  program. 

In  order  to  derive  closed-form  solutions  to  the  equa- 
tions of  motion  for  the  missile,  it  was  assumed  that  the  drag 
coefficient  was  constant  over  the  prediction  interval. 

This  assumption  will  introduce  errors  into  the  predictions 
However,  this  error  should  be  reduced  the  closer  the  missile 
approaches  to  its  terminal  point.   Near  the  terminal  point  of 
the  flight  the  greatest  error  will  be  in  the  predicted  time 
until  impact  since  the  Z  acceleration  will  be  the  greatest 
at  that  time.   This  should  not  greatly  affect  the  X  and  Y 
prediction  since  the  X  and  Y  accelerations  will  tend  to  be 
very  small  at  that  point. 

The  method  used  to  solve  the  equations  is  to  first  solve 
the  Z  motion  equation  for  the  time  at  which  the  missile  Z  is 
equal  to  the  target  Z.   The  value  of  time  thus  found,  which 
is  the  predicted  remaining  flight  time  of  the  missile,  is 
then  used  with  the  current  X  and  Y  positions  and  velocities 
to  compute  the  predicted  impact  point. 
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A.   METHOD  I  -  DRAG  PROPORTIONAL  TO  VELOCITY 

In  method  I  the  drag  force  on  the  missile  is  considered 
to  be  proportional  to  the  velocity  of  the  missile.   The  equa- 
tions for  this  case  are  linear  differential  equations. 

The  equations  of  motion  for  method  I  are 

mX  =  -K  X  (D-l) 

mY  =  -K2Y  (D-2) 


mZ  =  -K3Z  -  mG  (D-3) 


where 


K,  =  p  X  A  (D-4) 

1     '  '  r 


K2  =  p | Y | Ar  (D-5) 


K,  =  p  Z  A  (D-6) 

3     r   I   I   r 


•        •  • 

in  which  the  X,  Y,  and  Z  are  the  current  values  of  the 

velocities  of  the  missile,  p  is  the  atmospheric  density,  and 

2 

A   is  a  reference  area  of  .114  meter  . 

The  solution  for  the  Z  motion  equation  D-3  has  the  form 


where 


0  =  -Ae~at  -  Bt  +  C  (D-7) 


mZ  (0)     2r 

a  =  — ™ —  +  ni_G  (D_8) 

K3      K32 

B  =  p-  (D-9) 

K3 

C  =  A  +  Z  (0)  -  Z.  (D-10) 

m       t 

K3 

a  =  -i-     .  (D-ll) 

m 
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The  problem  now  becomes  the  solution  of  the  time-domain 
equation  of  the  Z  motion  for  the  time  at  which  the  missile  Z 
is  equal  to  the  target  Z.   The  C  coefficient  in  equation  D-10 
is  set  to  reflect  this  desired  equality  already. 

The  method  chosen  to  solve  this  equation  is  the  Newton- 
Raphson  method.   Upon  commencement  of  the  iterations  all  of 
the  required  coefficients  were  computed  using  the  current 
values  of  the  positions  and  velocities.   The  termination  test 
constant  used  to  terminate  the  iterative  process  was  0.1  sec- 
ond.  Whenever  one  value  of  time  had  changed  from  the  previous 
value  of  time  by  an  amount  less  than  or  equal  to  the  test 
increment,  the  iterative  process  was  terminated.   The  value 
of  0.1  seconds  was  chosen  for  the  test  value  since  regardless 
of  the  sample  interval  given  to  the  simulation  program  the 
missile  model  was  advanced  in  increments  of  0.1  seconds. 
Thus,  whenever  the  time  was  less  than  or  equal  to  the  minimum 
increment  that  the  missile  could  be  advanced,  the  iteration 
was  stopped. 

The  value  of  time  which  was  returned  from  the  Newton- 
Raphson  iteration  was  then  substituted  into  the  equations  for 
X  and  Y  impact  point  prediction.   These  equations  are 

m^™(0)  ft* 

X   =  X  (0)  +  — £ (1.0-e  Pr)  (D-12) 


p     m        K 


1 


where 


mY  ( 0 ) 

Y   =  Y  (0)  +  — £ (1.0-e"Yr)  (D-13) 

p    m        K-j 


Kl 

e  =  -±  (d-14) 

m 
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K 
Y  -  £  .  (D-15) 

These  are  the  equations  for  method  I  prediction.   The 
method  II  equations  are  more  complicated  to  solve  since  they 
are  non-linear. 

B.   METHOD  II  -  DRAG  PROPORTIONAL  TO  VELOCITY  SQUARED 

In  this  derivation  of  impact  prediction  equations  the 
drag  force  on  the  missile  is  taken  to  be  proportional  to  the 
square  of  the  missile  velocity.   Again,  the  drag  coefficients 
are  computed  at  the  beginning  of  the  solution  and  considered 
constant  throughout. 

The  equations  of  motion  are 

mX  =  -I^X2  (D-16)  ■ 

mY  =  -K2Y  (D-17) 

mZ  =  -K3Z2  -  mG   .  (D-18) 

The  steps  of  this  solution  are  the  same  as  used  in  the 

previous  method:   1)  solve  the  Z  motion  equation  for  Z  as  a 

function  of  time,  2)  solve  for  the  time  when  missile  and 

target  Z  are  equal,  and  3)  use  the  time  thus  found  to  compute 

X   and  Y   from  the  prediction  equations. 
P      P  F  * 

First  solve  for  Z  as  a  function  of  time,  writing  equation 
D-18  as 

|£  =  -(Cy2+G)  (D-19) 

where 

K3 
C  =  -=■  (D-20) 

m 

and 

■ 

y  =  Z    .  (D-21) 
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Separating  variables  and  integrating, 

1      -1 


/CG 


tan   (y/C/G)  +  D,  =  t  .  (D-22) 


At   t   =    0,    y   =    Z     (0).       Therefore 

D,    =    -=-  tan_1(Z     (0)/c7G)        .  (D-23) 

1         /CG 

Solving    for   y    and   substituting    for   D, , 

y    -    /G7C    tan^/CG)    =    -/G/C    tan(t/CG)  (D-24) 

or 

• 

y    =    -/G7C    tan(t/CG)    +    Z     (0)        .  (D-25) 

m 

dZ 
Substituting   y   =   3—  and   letting    Z    (0)    be    constant. 

3   J         dt  3      m 


and 


dZ    =     [-/G7C    tan(t/CG)    +    Z     (0)]dt  (D-26) 

m 


1 


Z    =    -/G/C    -=-   log      cos(t/CG)    +    Z     (0)t    4-    D0     .  (D-27) 

/CG  e  m  2 

At    t   =    0 ,    Z    =    Z     (0) .      Therefore    Dn    =    Z     (0)     and 

m  2  m 

Z(t)    =    -i-  log      cos(t/CG)    +    Z     (0)t    +    Z     (0)  (D-28) 

L  e  mm 

where    Z    (0)    and   Z    (0)    refer   to   the    Z    velocity   and  position 

mm  j.  - 

of  the  missile  at  the  time  the  prediction  is  made. 

Solving  for  the  time  when  Z (t)  is  equal  to  the  target 
elevation  is  done  by  substituting  Z (t)  =  Z  and  using  the 
Newton-Raphson  method. 

The  solution  of  the  X  and  Y  motion  equations  for  impact 
point  will  be  shown  for  the  X  equation  only  and  the  Y  equa- 
tion will  be  written  by  inspection. 

mX  =  -K^2  (D-29) 


&  =  -Cy2  (D-30) 


dt 
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K 

where  y  =  X  and  C  =  — .   Separating  variables  and  integrating 


m 
1 


+  D1  =  t  .  (D-31) 

At  t  =  0.  y  =  X  (0) .   Therefore 

D   = s-i .  (D-32) 

cxm(o) 

dX 
Substituting  y   =    -?r-  and   letting  X    (0)    be    constant 
3    -*         dt  3      m 

§$•  =    Ct    -   CD,       .  (D-33) 

dX  1 

Separating  variables  and  integrating 

X  =  i-  log  (Ct  -  -r-^ — )  +  D9   .  (D-34) 

XtJ0> 
m 

At  t=  0,  X=  X  (0).   Therefore 

D,  =  X  (0)  -  i  log  (-  -^ )  (D-35) 

2     m      C    ye      (Q) 

m 


and 


X(t)  -'  i  log  (Ct  -  -r-i — )  -  £  log  (-  — ^ — )  +  X  (0)  (D-36) 

L    e      X  (0)     L    e    X(0)  m 

m  m 

X(t)  =  h   log  (Ct  -  T-i )  (-X(0))  +  Xm(0)      (D-37) 

C    e      X  (0)     m        m 

m 


or 

K,tX  (0) 

X(t)  =  g-  log  [1  -  -^— £ ]  +  Xm(0)  .       (D-38) 

K,    3e         m        m 

For  this  method  the  prediction  equations  are 

K,tX  (0) 

X   =  §-  logo[l  -  -=—. 2 ]  +  X(0)         (D-39) 

p   K,    3e         m        m 


and 


K9tY(0) 

Y^  =  £-  log  [1  -  -2— £ ]  +  Y(0)   .     (D-40) 

p   K~    3e         m        m 
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APPENDIX  E 
REFERENCE  SYSTEMS  FOR  THE  PROBLEM 

The  reference  systems  for  the  problem  are  based  on  the 
position  of  the  launch  point  of  the  aircraft  involved.   The 
inertial  navigator  of  the  aircraft  is  initialized  when  the 
aircraft  is  launched.   All  other  points  in  the  problem  are 
also  referred  to  this  initialization  point  of  the  aircraft 
even  though  they  may  not  be  measured  with  respect  to  that 
point  at  first.   The  measured  missile  state  vector  is  an 
example  of  this  type  in  that  it  is  first  referred  to  the 
mother  aircraft  and  then  converted  to  a  vector  which  is 
referred  to  the  initialization  point. 

The  cartesian  reference  is  a  standard  three-dimensional, 
right-hand  reference  system.   Whenever  a  reference  is  shifted, 
the  new  frame  is  also  a  three-dimensional,  right-hand  system. 
Figure  22  illustrates  the  various  reference  points  and  systems 

Azimuth  angle  is  defined  as  the  angle  between  true  north, 
positive  clockwise  to  the  line-of-sight  from  the  aircraft  to 
the  item  of  interest.   In  order  to  define  the  depression 
angle,  it  is  first  necessary  to  define  the  normal  plane.   The 
normal  plane  is  defined  as  being  a  plane  normal  to  a  line  pas- 
sing from  the  aircraft  to  the  center  of  the  earth.   Depression 
angle  is  defined  as  the  angle  between  the  normal  plane  and  the 
line-of-sight.   This  LOS  must  be  contained  in  a  plane  which  is 
perpendicular  to  the  normal  plane.   These  angles  are  shown  in 
Fig.  23  and  24  respectively. 
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Figure    22.      Problem  Cartesian   Reference    Frame 
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APPENDIX  F 


THE  MISSILE  MODEL 


The  missile  model  was  relatively  honest  in  its  aerodynamic 
performance  although  it  was  not  a  true  missile.   The  maximum 
yaw  rate  of  the  missile  was  limited  to  allow  radial  accelera- 
tion loads  not  to  exceed  ten  G's.   This  was  thought  to  be 
representative  of  current  capabilities. 

The  missile  was  powered  by  a  rocket  motor  of  3000-lbs. 
thrust  with  a  burn  time  of  30  seconds.   After  the  burnout  the 
missile  glided  for  the  remainder  of  its  flight. 

The  model  was  so  designed  that  the  mass  and  moment  of 

inertia  changed  as  the  rocket  motor  burned.   The  total  pro- 

pellant  weight  was  approximately  163.3  kg.   The  reduced  mass 

and  moment  of  inertia  were  then  used  for  the  remainder  of  the 

flight.   A  free  body  diagram  of  the  missile  appears  in  Fig. 

25.   The  equations  of  motion  are 

mR  =  -D  cosy  -  L  siny  +  T  cos6  (F-l) 

mZ  =  -D  siny  +  L  cosy  +  T  sine  -  mG       (F-2) 

I   9  =  -jL  -  jD  sinct  +  ICF    -    K.Q0         (F-3) 
yy       J      J  4 


where 


Q  =  pV2  (F-4) 

D  =  1^0  (F-5) 

L  =  KcQa  (F-6) 
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and 
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K. 

Kt 
K, 


m/sec 


thrust) 

lift) 

drag) 

velocity) 

control  force) 

velocity  vector  angle)  rad 

angle  of  attack)        rad 

missile  axis  angle) 

axis  rate) 

axis  acceleration) 

mass) 


13344.  N 

N 


N 


N 


gravity) 

control  moment  arm) 

center  of  pressure 
moment  arm) 

moment  of  inertia) 

air  density) 
drag  constant) 

control  force 
constant) 

damping  constant) 
lift  constant) 


rad 
rad/sec 

rad/sec' 


4  5  3.6  kg 

(290.3  kg  after  burnout) 

9.8  m/sec2 

2  .  m 

1.  m 

34  0.2  nt-m2 

2 
(217.7  nt-m   after  burnout) 

1.2  84  kg/m3 

.00556  nt/ (kg/m3) (m/sec)2 


3         2 

.01  nt/(kg/m  ) (m/sec)  rad 

.004  nt/ (kg/m3) (m/sec) 2 (rad/sec) 

3         2 
.008  nt/(kg/m)  (m/sec)  rad. 


71 


APPENDIX  G 
THE  LAUNCHING  AIRCRAFT  MODEL 

The  launching  aircraft  model  was  a  very  simple  model  for 
the  simulation.   The  aircraft  was  given  an  initial  heading 
and  velocity.   It  was  then  assumed  to  fly  a  straight  line 
making  no  turns.   The  location  of  the  aircraft  was  updated 
each  sample  period. 

Although  preliminary  investigation  had  shown  that  the 
most  critical  component  of  noise  in  the  entire  system  was 
that  in  the  inertial  navigation  package,  the  inertial  navi- 
gator was  considered  to  be  noiseless.   The  reasons  for 
neglecting  the  noise  were  that  1)  the  primary  purpose  of  this 
investigation  was  in  the  field  of  application  of  impact  pre- 
diction to  an  operator-controlled,  air-to-surface  weapon 
system,  and  2)  no  figures  were  available  concerning  sample- 
to-sample  noise  in  inertial  navigators. 
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APPENDIX  H 
MAN/SYSTEM  INTERFACE 

In  any  system  with  a  human  operator,  it  is  important  that 
information  be  presented  concisely  and  in  such  a  manner  that 
it  can  best  be  used  by  the  operator.   The  communication  link 
in  any  such  application  is  a  two-way  channel.   The  system 
sends  information  to  the  operator  in  one   direction  and  the 
operator-generated  control  commands  travel  in  the  other. 

The  link  to  the  operator  carried  two  items  of  information 
The  first  was  the  most  recent  prediction  of  impact  point  and 
the  second  was  the  predicted  time  until  impact. 

The  first  item,  impact  point  prediction,  was  displayed 
to  the  operator  on  a  display  scope  with  the  target  assumed  to 
be  in  the  center.   This  presentation  was  limited  in  such  a 
manner  that  the  predicted  impact  point  could  never  appear 
off-scale.   Figures  26,  27,  and  28  give  examples  of  the 
limiting  action  which  was  employed.   The  display,  when  in  a 
limited  state,  gave  no  information  as  to  the  error  of  the 
predicted  impact  point  except  the  direction  in  which  it  was 
off-scale. 

The  display  was  generated  such  that  the  missile  always 
approached  the  target  from  the  bottom  of  the  scope.   This 
required  a  transformation  of  the  axes  of  the  display.   This 
was  necessary  so  that  the  motion  of  the  control  stick  always 
drove  the  missile  toward  the  target. 
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Figure  26.   Impact  Point  Within  Scale 
on  Display 


Figure  27.   Impact  Point  Off-Scale 
Left  on  Display 


Figure  28. 


Impact  Point  Off-Scale 
Left,  Off-Scale  Vertically 
on  Display 
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The  transformation  of  axes  consisted  of  both  a  translation 
and  rotation.   The  translation  was  accomplished  by  generating 
a  reference  frame  which  had  its  origin  located  on  the  target. 
The  predicted  impact  point  was  then  referred  to  this  origin 
and  the  axes  rotated  in  such  a  manner  as  to  make  the  Y  axis 
of  the  rotated  system  coincide  with  the  LOS  from  the  current 
measured  missile  position  to  the  target.   This  process  is 
illustrated  in  Fig.  29  and  30. 

It  can  be  shown  that  the  location  of  the  predicted  im- 
pact point  in  the  translated  and  rotated  reference  frame 
(display  reference  system)  is 


X,  =  AX,- 

— d     — df 


(H-l) 


where 


*d  = 


X 


(H-2) 


A  = 


cos  6   -sin  0 


sin  0    cos  0 


5df 


X  -  X, 

p   t 

Y  -  Y, 

P     t 


(H-3) 


(H-4) 


and 


0  =  tan_1((X.  -X J/(Y,  -Y  J  )   . 
t  m    t  m 


(H-5) 


The  second  item  of  information  presented  to  the  operator 
was  the  predicted  time  remaining  until  impact.   This  value 
was  computed  once  each  sample  period  throughout  the  missile 
flight.   It  was  generated  on  the  digital  computer,  converted 
to  an  analog  voltage,  and  displayed  on  a  digital  voltmeter. 
The  display  was  scaled  such  that  times  were  displayed  only  if 
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+  X' 


^MiM 


Figure  29.   Translation  of  Axes  for  Display 


+Y' 


I^MM 


Figure  30.   Rotation  of  Axes  for  Display 
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they  were  less  than  or  equal  to  100  seconds.   Any  value  of 
time  which  was  greater  than  100  seconds  was  displayed  as  100 
seconds . 

Commands  from  the  operator  to  the  system  formed  the  other 
direction  of  the  two-way  communication  link.   These  commands 
were  generated  by  a  control  stick  which  the  operator  used  to 
guide  the  missile.   Full  deflections  on  the  stick  corresponded 
to  full  deflections  on  the  missile  controls.   The  pitch  con- 
trol was  scaled  to  provide  45  degrees  of  control  movement  on 
the  missile  when  the  stick  was  fully  deflected.   The  yaw  con- 
trol provided  a  yaw  rate  of  20  degrees/second  when  fully  de- 
flected.  This  value  was  also  subject  to  limiting  by  the 
missile  dynamics. 
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APPENDIX  I 
A  DISCUSSION  OF  THE  SIMULATION  PROGRAM 

Tnis  appendix  describing  the  program  used  in  the  simula- 
tion of  this  problem  is  included  in  the  event  that  some  other 
investigator  would  desire  to  continue  research  in  this  oar- 
ticular  field.   In  it  the  author  has  tried  to  include  all 
pertinent  information  concerning  the  routines  and  their  func- 
tions and  especially  the  programming  "tricks"  which  were  used 
in  some  instances  to  provide  a  particular  function  in  the 
overall  program. 

The  appendix  is  divided  into  four  sections.   The  first 
section  is  a  summary  of  the  routines  and  the  function  which 
each  performs.   The  second  section  lists  the  routines  and 
the  variable  names  within  each  routine.   The  method  of  inter- 
routine  communication  is  the  subject  of  the  third  section 
and  the  programming  "tricks"  used  are  discussed  in  a  fourth 
section.   A  complete  listing  of  the  program  and  an  analog 
patching  diagram  follow  the  appendix. 

A.   SUMMARY  OF  THE  ROUTINES  AND  THEIR  FUNCTIONS 
1 .   Main  Line  Program 

The  main  line  program  serves  to  perform  two  functions. 
The  first  is  as  a  vehicle  for  input  to  the  program  and  output 
from  the  program,  and  the  second  is  as  a  calling  and  control- 
ling sequence  for  the  various  routines  which  do  the  actual 
computation  for  the  simulation.   If  operated  in  a  segmented 
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mode  of  core  allocation ,  the  main  line  program  is  always  core 
resident. 

Several  test  functions  are  performed  by  the  main  line 
program  also.   It  tests  for  the  condition  of  problem  termi- 
nation.  Problem  termination  occurs  when  the  measured  Z  of 
tne  missile  is  equal  to  or  less  than  the  target  Z.   It  also 
controls  switching  from  the  passive  target  location  mode  of 
the  problem  to  the  active  mode  in  which  the  missile  is  fired 
and  guided  to  the  target. 

2 .   Initialization  Routines 

a.  Inertial  Navigator  Initialization 

The  subroutine  INIT  initializes  the  inertial 
navigation  routine  by  setting  the  initial  position  of  the 
aircraft  including  the  altitude.   It  also  initializes  the 
heading  and  computes  component  velocities .   The  final  items 
initialized  by  the  routine  are  the  elapsed  time  counter, 
missile  flight  time  counter,  and  the  problem  time  increment. 

b.  Missile  Initialization 

The  subroutine  MINIT  initializes  the  missile  model 
at  the  time  of  launch.   All  parameters  of  the  missile  model 
are  set  equal  to  the  aircraft  parameters  at  the  time  of  launch 
Also  the  routine  gives  initial  values  to  the  measured  missile 
vector  to  enable  other  routines  to  function  properly  during 
the  first  second  of  the  simulation  until  all  routines  are 
generating  full  outputs. 
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3 .  Aircraft  Simulation  Routine 

The  subroutine  NAVPCK  is  the  routine  which  provides 
aircraft  state  vector  information  for  the  problem.   The 
routine  provides  a  straight-line  flight  path  for  the  aircraft 
at  any  initial  heading  from  0  to  359  degrees  at  any  initial 
velocity.   The  routine  receives  flight  path  information  in 
R/0  form  and  provides  flight-path  information  in  the  three 
cartesian  coordinates. 

4 .  Passive  Target  Location 

The  subroutine  PFIL  simulates  the  action  of  the  hybrid 
filter  presented  in  section  II-A-3.   This  location  action  is 
simulated  by  interpolating  and  normalizing  values  of  variance 
provided  for  the  hybrid  method  and  generating  noise  which  is 
added  to  the  true  target  location.   The  final  noisy  location 
of  the  target  is  then  used  throughout  the  remainder  of  the 
problem. 

5 .  Missile  Model 

The  missile  model  routine  is  of  sufficient  importance 
that  an  entire  previous  appendix  was  devoted  to  it.   The  reader 
is  referred  to  Appendix  F  for  a  discussion  of  the  missile  model- 

6 .  Routines  Involved  in  Missile  Tracking 
a.   Radar  Vector  Generation 

The  subroutine  REVCON  is  the  routine  which  gener- 
ates the  true  radar  vector  from  the  geometrical  relationships 
between  the  true  missile  and  the  aircraft.   This  vector  reo- 
resents  the  true  radar  state  before  any  noise  is  added  and 
before  quantization  takes  place. 
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b.  Noise  Addition  and  Quantization  of  Radar  Vector 
The  subroutine  RADAR2  takes  the  true  radar  vector 

produced  by  REVCON  and  performs  two  operations  on  it.   First 
the  routine  adds  the  measurement  noise  associated  with  the 
particular  measured  state  variable.   Second,  the  routine 
quantizes  the  measurements  into  the  discrete  form  which 
would  be  provided  by  a  digitized  radar  output  unit. 

c.  Radar  Vector  Filtering 

The  subroutine  AFIL  provides  the  filtering  of 
tiie  noisy,  quantized  radar  vector  to  produce  a  radar  vector 
whose  variance  is  less  than  that  of  the  noisy  vector.   The 
subroutine  uses  a  schedule  of  precomputed  filter  gains  and 
a  switch  which  determines  the  number  of  samples  after  which 
the  gain  schedule  is  truncated.   This  value  can  be  from  one 
to  thirty  seconds  assuming  gains  computed  for  a  sample  time 
of  one  second. 

d.  Spherical  to  Cartesian  Coordinate  Conversion 
The  subroutine  CORDCN  processes  the  filtered 

radar  vector  spherical  values  to  produce  a  measured  missile 
vector  composed  of  missile  position  and  velocities  relative 
to  the  true  problem  reference.   The  subroutine  first  pro- 
duces a  cartesian  vector  relative  to  the  aircraft.   The  air- 
craft component  positions  and  velocities  are  then  added  to 
the  relative  components  to  produce  the  final  measured  missile 
vector. 
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7 .  Initial  Missile  Control  Routine 

The  subroutine  SHOOT  controls  the  flight  of  the  mis- 
sile during  the  boosted  phase  (first  thirty  seconds)  of  its 
flight.   The  routine  attempts  to  drive  the  missile  velocity 
vector  to  assume  an  up  angle  of  4  7.5  degrees  and  to  turn  the 
missile  velocity  vector  projection  on  the  XY  plane  to  coin- 
cide with  the  LOS  from  the  missile  to  the  target. 

8 .  Impact  Prediction  and  Display 

a.  Impact  Point  Prediction 

The  subroutine  PREDCT  develops  the  X  and  Y  im- 
pact point  prediction  values.   This  routine  employs  the  sub- 
routine NR  to  assist  in  this  task. 

b.  Time-Until-Impact  Solution 

The    subroutine   NR   solves    for   the    time    to   go 
until    impact    for   the   missile   by    using   the   Newton-Raphson 
iterative   method  of   solution    for   differential    equations. 

c.  Command  Generation 

The  subroutine  COMGEN  converts  the  outputs  of 
the  control  stick  to  values  which  can  drive  the  missile  mod- 
el.  This  routine  uses  the  SDS  9300  routines  for  analog-to- 
digital  conversion.   Full-scale  values  of  pitch  correspond 
to  plus  or  minus  45  degrees  of  control  deflection  and  full- 
scale  values  of  yaw  correspond  to  turn  rates  of  plus  or 
minus  20  degrees/second. 

9 .  Prediction  Display  and  Statistics  Routines 
a.   Display  Generation 

The  subroutine  DISGEN  converts  the  predicted  im- 
pact points  into  values  which  can  be  used  to  drive  the 
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situation  display.   This  routine  performs  the  axis  rotation 
and  translation  necessary  to  generate  a  useful  display.   The 
limiting  of  the  display  range  is  also  performed, 
b.   Minimal  Statistics  Routine 

The  subroutine  MISS  provides  minimal  statistics 
about  the  success  of  the  run  at  the  time  of  flight  termina- 
tion.  This  routine  provides  the  X  and  Y  miss  distances, 
impact  time  of  the  missile,  and  total  miss  distance.   More 
complete  statistics  are  generated  by  a  separate  program  which 
analyzes  data  which  is  written  on  a  storage  medium  during  the 
execution  of  the  program.   The  separate  analysis  routine  pro- 
vides such  results  as  missile  flight  profile,  control  applied 
vs.  time,  missile  location  errors  vs.  time,  and  others. 


B.   ROUTINE  SYMBOLS  AND  THEIR  MEANINGS 


1 .   Main  Line  Program 
Symbol 
GBQ/KBQ 

RUNS/RNS 

NP 

SMU,    ARNG,    ARNGD,    AAZ ,    ADP 


VL,    THI,     Z,    TME 


ETIME,  TIME,  FTIME 


ETA 


Meaning 

Time  for  active  filter 
gain  truncation. 

Number  of  runs  for  the 
program. 

Switch  used  to  initialize 
AFIL. 

Mean  of  noise  and  noise 
variances . 

Initial  velocity,  heading 
(in  degrees) ,  aircraft 
altitude,  and  time  incre- 
ment. 

Elapsed  time,  time  incre- 
ment, flight  time. 

Initial  launch  angle. 
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Subroutine  INIT 


Symbol 

XA , Y A , Z  A , XDA , YDA , Z  DA 

CK,CT,TH,V 


3.   Subroutine  NAVPCK 
Symbol 
CK,CT 


XA , Y A , Z  A , XDA , YDA , Z  D A , TH , V 

4 .  Subroutine  MINIT 
Symbol 

YAW, PITCH 

XMM , YMM ,  ZMM , XDMM , YDMM , 
ZDMM 

XM , YM , Z  M , XDM , YDM , Z  DM , 
XDDM,YDDM,ZDDM, 
HEAD , GAMMA , THET , 
THETD,THETDD 

5 .  Subroutine  PFIL 
Symbol 

X,Y 
QKfQM^^J 

XS,YS,A 
XT,YT,CA 


Meaning 

Aircraft  location  and 
velocity  in  cartesian 
coordinate  frame. 

Conversion  factor  (degrees 
to  radians) ,  converted 
heading  angle,  heading  in 
degrees ,  velocity  in 
meters/ second. 


Meaning 

Conversion  factor  (degrees 
to  radians) ,  converted 
heading. 

Aircraft  state  vector. 


Meaning 

Missile  command  vector. 

Measured  missile  vector 
in  position  and  velocity 

True  missile  vector. 


Meaning 

True  target  location. 

Constants  used  for 
interpolation. 

X  variance,  Y  variance, 
noise  value  returned  by 
GAUSS. 

Noisy  target  location, 
CPA  in  problem. 
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6.   Subroutine  MISAD 


Symbol 
RHO,T,G 

RE  FAR 

K1,K3,K4 ,K5 

K,YINC,ANG,VEL,VELSQ 


DELTA,  GAMMA,  ALPHA 

D,QL,QM 

QMASS , QI YY , QMULT 

7.   Subroutine  REVCQN 
Symbol 
X,Y,Z 
RPRI ,RDPRI ,T1 ,T2 ,VLXY 


R,RD,T,P 

8. 

Subroutine  RADAR2 

Symbol 

CI 

I TEMP, AN 

R , RD , T  ,  P 


Meaning 

Air  density,  thrust, 
gravity. 

Aerodynamic  reference 
area. 

Drag,  lift  and  moment 
constants . 

Advance  count,  yaw  incre- 
ment, current  missile 
heading,  current  velocity, 
current  velocity  squared. 

Control  angle,  missile 
axis  angle,  angle  of 
attack. 

Drag,  lift,  and  moment 
forces . 

Missile  mass ,  moment  of 
inertia,  constant  relating 
mass  and  moment. 


Meaning 

Work  variables. 

Work  variables  for  develop- 
ing the  radar  vector. 

True  radar  vector. 


Meaning 

Bin  size  for  13-bit  shaft 
encoder. 

Temporary  variable  used 
in  the  conversions,  and 
noise  value  returned  by 
GAUSS. 

Noisy,  quantized  radar 
vector. 
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9 .       Subroutine   AFIL 
Symbol 
NP 

ZR1,ZR2,ZA1,ZD1 

XR1 , XR2 , XA1 , XA2 , XD1 , XD2 


GRG,GAZ,GDP 

10.   Subroutine  CQRDCN 
Symbol 
CSP,SNP,CST,SNT 


CSPSNT , CSPCST , SNPSNT , 
SNPCST 

R,RD,T,TD,P,PD 


XMM,YMM,ZMM 

XDMM,YDMM,ZDMM 

11.   Subroutine  SHOOT 
Symbol 
GAMMA,  HEAD 

ETA,A1 


YAW, PITCH 


Meaning 

Switch  to  control  filter 
initialization. 

Measurements  of  range, 
range  rate,  azimuth,  and 
depression  angles. 

Estimated  values  of  range, 
range  rate,  azimuth, 
azimuth  rate,  depression, 
and  depression  rate. 

Range , azimuth ,  and  de- 
pression filter  gains. 


Meaning 

Sine  and  cosine  of  (f> 
and  8 . 

Cross  products  of  sine 
and  cosine  values. 

Range,  ranqe  rate,  theta, 
theta  rate,  phi,  and  phi 
rate. 

Measured  missile  cartesian 
position. 

Measured  missile  cartesian 
velocity. 


Meaning 

Current  missile  velocity 
vector  elevation  and  cur- 
rent missile  heading. 

Desired  missile  velocity 
vector  elevation  and  cur- 
rent LOS  angle  from  the 
missile  to  the  target. 

Missile  command  vector. 


86 


12.   Subroutine  PREDCT 


Symbol 

QMASS ,RHO,REFAR 

K2  ,K3,V 
XP,YP,ZP 

13.   Subroutine  NR 
Symbol 
QMASS , RHO , REFAR ,  G ,  V 


Kl, ALPHA 
A,B,C,TOLD,TNEW 

FT , FDT 

14.  Subroutine  COMGEN 
Symbol 

YAW, PITCH 
Numeric  arguments 

15.  Subroutine  DISGEN 
Symbol 

XDF,YDF 

XTA,YTA 
ALPH 


Meaning 

Missile  mass,  air  density/ 
missile  aerodynamic  refer- 
ence area. 

Drag  coefficients  and 
missile  velocity. 

X  and  Y  impact  position 
predictions,  and  time 
until  impact. 


Meaning 

Missile  mass,  air  density, 
missile  aerodynamic  refer- 
ence area,  gravity,  and 
missile  velocity. 

Drag  coefficient  and 
exponential  argument. 

Coefficients  of  prediction 
equations ,  last  time  esti- 
mate ,  new  time  estimate. 

Function  value ,  and 
derivative. 


Meaning 

Missile  Command  vector 

Conversion  factors. 

Meaning 

Translated  predicted 
impact  point. 

Used  to  determine 
rotation  angle. 

Rotation  angle. 
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XD,YD,ZD 

16.   Subroutine  MISS 
Symbol 
ZTM1,TI1,TI2 

XIM1,YIM1,XI,YI 


ETI,FTI 
Kl,K2,QK2 

XMS , YMS , XT , YT 


Display  points  and  time 
until  impact. 


Meaning 

Last  preimpact  missile  Z 
position,  times  which 
bracket  impact  time. 

Last  preimpact  X  and  Y 
positions,  and  missile 
impact  positions  in  X 
and  Y. 

Impact  time  as  elaosed 
time  and  flight  time. 

Variables  used  to  edit 
inaccuracies  out  of 
flight  time  clock. 

X  and  Y  miss  distances, 
and  last  estimated  target 
position. 


C.   INTERROUTINE  COMMUNICATION 

Interroutine  communication  is  carried  out  throughout  th« 
program  by  means  of  labeled  FORTRAN  COMMON  blocks.   These 
blocks  are  designed  in  such  a  manner  that  variables  relatincr 
closely  to  each  other  have  been  grouped  into  sections  named 
to  reflect  their  functions.   The  block  names  are  given  below 
along  with  the  variables  included  and  a  description  of  the 
function  of  the  block. 

1.   Blank  COMMON 

The  variables  included  in  blank  COMMON  are  IX,  ETIME, 
TIME,  and  FTIME.   IX  was  originally  the  seed  value  for  the 
GAUSS  routine  which  produces  random  numbers .   It  was  not  re- 
tained when  a  change  in  random  number  generators  was  made. 
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ETIME,  TIME,  and  FTIME  are  the  problem  elapsed  time,  time 
increment,  and  missile  flight  time  respectively. 

2.  Block  TGPRD 

The  variables  included  in  the  block  TGPRD  are  the 
cartesian  location  of  the  target  and  the  X  and  Y  predicted 
impact  points  of  the  missile.   ZP  is  the  time  until  impact. 

3.  Block  SIGMAS 

This  block  contains  the  values  of  the  sigmas  to  be 
used  in  the  radar  to  introduce  noise  in  the  measurements. 
These  values  are  set  in  the  main  line  program. 

4.  Block  GAINS 

The  computed  filter  gains  used  by  the  active  radar 
filter  are  stored  in  this  dimensioned  block.   The  block 
provides  storage  for  120  range  filter  gains  and  60  azimuth 
and  depression  filter  gains  each. 

5.  Block  MISLM 

This  block  contains  the  values  of  measured  missile 
position  and  velocity  developed  by  the  missile  tracking 
routine . 

6.  Block  ACRFT 

The  values  of  aircraft  position  and  velocity  at  each 
time  k  are  held  in  this  block.   Also  this  block  contains  the 
aircraft  velocity  in  meters  per  second  and  aircraft  headinq 
in  degrees . 

7.  Block  MISLT 

This  vector  contains  all  values  which  pertain  to  the 
true  missile.   The  vector  contains  position,  velocity,  and 
acceleration  in  cartesian  coordinates;  missile  axis  angle, 
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angular  velocity,  and  angular  acceleration;  and  velocity  vec- 
tor heading  and  elevation  angle.  This  is  updated  each  sample 
period  to  reflect  the  true  state  of  the  missile. 

8.  Block  COMND 

The  commands  to  the  missile  in  yaw  and  pitch  comprise 
the  two  elements  of  this  block.   The  yaw  values  are  up  to 
twenty  degrees  per  second  in  either  direction  and  pitch  val- 
ues range  from  plus  to  minus  forty -five  degrees  in  control 
deflection. 

9.  Block  SIGS 

The  values  of  variance  which  are  used  in  the  passive 
location  routine  PFIL  are  transferred  to  the  routine  by  this 
block.   Each  of  the  values  is  dimensioned  to  accept  fifty- 
nine  values  of  variance. 
10.   Block  RADAR 

This  block  is  used  to  transfer  the  values  of  the  radar 
information  through  all  of  the  routines  of  the  missile  track- 
ing portion  of  the  program.   The  six  elements  of  the  vector 
are  range,  range  rate,  azimuth,  azimuth  rate,  depression,  and 
depression  rate. 

D.   PROGRAMMING  "TRICKS" 

Only  one  programming  irregularity  exists  in  the  simula- 
tion program.   In  the  subroutine  PFIL  the  variable  ZP  in  the 
COMMON  block  TGPRD  is  used  as  a  switch  to  allow  the  routine 
to  hold  the  true  target  location  after  the  first  call.   The 
variable  ZP  was  chosen  since  it  was  already  in  the  COMMON 
block  and  since  it  is  not  used  by  any  other  routine  in  the 
passive  phase  of  the  program. 
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